**FULL TITLE** 

ASP Conference Series, Vol. **VOLUME**, **YEAR OF PUBLICATION** 
**NAMES OF EDITORS** 



Keeping the Universe ionised: photoheating and the 
high-redshift clumping factor of the intergalactic gas 



Andreas H. Pawlik, 1,2 Joop Schaye, 2 and Eveline van Scherpenzeel 2 

Abstract. The critical star formation rate density required to keep the in- 
tergalactic hydrogen ionised depends crucially on the average rate of recom- 
binations in the intergalactic medium (IGM). This rate is proportional to the 
clumping factor Cigm = (pt)iGM/ 'p\, where pb and pb are the local and cosmic 
mean baryon density, respectively and the brackets ()igm indicate spatial aver- 
aging over the recombining gas in the IGM. We perform a suite of cosmological 
smoothed particle hydrodynamics simulations to calculate the volume-weighted 
clumping factor of the IGM at redshifts z > 6. We investigate the effect of 
photoionisation heating by a uniform ultraviolet background and find that pho- 
toheating strongly reduces the clumping factor as the increased pressure support 
smoothes out small-scale density fluctuations. Even our most conservative esti- 
mate for the clumping factor, Cigm = 6, is five times smaller than the clumping 
factor that is usually employed to determine the capacity of star-forming galaxies 
to keep the z « 6 IGM ionised. Our results imply that the observed population 
of star-forming galaxies at z w 6 may be sufficient to keep the IGM ionised, 
provided that the IGM was reheated at z ;> 9 and that the fraction of ionising 
photons that escape star-forming regions to ionise the IGM is larger than 0.25. 



1. Introduction 

The absence of a Gunn-Peterson trough in many of the absorption spectra to- 
wards high-redshift quasars suggests that the re ionisation of intergalactic hy - 
drogen was completed at redshifts z ^ 6 (e.g., iFan. Cariili.&; Keatingl [2006). 



Current observational estimates of the ultraviolet (UV) luminosity density a t 
redshifts z <J 6 (for a comprehensive discussion see, e.g., iBouwens et al.ll2007l ). 



on the other hand, may imply star formation rate (SFR) densities several times 
lower than the critical SFR density required to keep the intergalactic medium 
(IGM) ionised. 

The critical SFR density, 
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(1) 



originally derived bv lMadau. Haardt. k. Reesl ( 19991 ) and here re scaled to match 



the m ost recent WMAP estimate for the cosmic baryon density (jKomatsu et al 



2008j), is inversely proportional to the escape fraction / esc , i.e. the fraction of 
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ionising photons produced by star-forming galaxies that escape the interstellar 
medium (ISM) to ionise the IGM. It is proportional to the average recombination 
rate in the IGM, parametrised using the dimensionless clumping factor Cigm = 
(/°b)iGM/Pb) where pb is the baryon density, is the mean baryon density of 
the Universe and the brackets ()igm indicate spatial averaging over the gas 
constituting the recombining IGM. 

Most observational studies that compare the SFR density derived from esti- 
mates of the UV luminosity density at redshift z ~ 6 to the critical SFR density 
assume an escape fraction / esc <^ 0.5 and a clumping factor Cigm = 30. Con- 
sequently, the observed population of galaxies has been found to be incapable 
of keeping the intergalactic hydrogen ionised, forming massive stars at a rate 
which is up to an order of magnitude lower than required by Eq. [TJ It has, 
however, been noted that the discrepancy between the observationally inferred 
and critical SFR densities at z ~ 6 could be resolved if the employed clumping 
factor were too high (e.g., ISawicki Thompsonll2006l ). 

In this work we perform cosmological Smoothed Particle Hydrodynamics 
(SPH) simulations to compute the clumping factor of the IGM. We pay partic- 
ular attention to the fact that the clumping factor depends on the definition of 
which gas comprises the IGM. By comparing the critical SFR density, updated 
using our estimate of the clumping factor, with the observationally inferred 
SFR density we argue that the observed population of UV galaxies may well 
be capable of keeping the z ~ 6 Universe ionised. We also investigate the ef- 
fect of photoionisation heating by a uniform UV background on the evolution 
of the clumping factor. We demonstrate that photoheating significantly lowers 
the clumping factor. Our results are insensitive to the redshift at which the UV 

background is turned on. 

This work is described in more detail in lPawlik. Schave. van Scherpenzeel 

(120091 1. 



Simulations 



We u se a modified version of the N-body/TreePM/SPH code gadget-2 (jSpringe] 



2005) to perform a suite of cosmological SPH simulations including radiative 



cooling. We assume a flat ACDM universe and employ the set of cosmological 
parameters [U m , O bj Q A , cr 8 , n s , h] given by [0.258 0.0441, 0.742 0.796 , 0.963, 
0.719], consistent with the WMAP 5-year results (jKomatsu et al.ll2008l ). 



The gas is of primordial composition, with a hydrogen mass fraction X = 
0.752 and a helium mass fraction Y = 1 — X. Radiative cooling and heating are 
i nclud ed assuming ionisation equilibrium, as described in lWiersma. Schave. &: Smith 



(2009). Molecular hydrogen and deuterium and their catalysts are kept photo- 



disso ciated by a soft UV backgrou nd and never contribute to the cooling rate 
(e.g., lHaiman. Rees. fe Loebl fl997l ) . Star form ation is modelled by turning g as 
particles into star particles using the recipe of ISchave fc Dalla Vecchial (|2008l ). 



All simulations use 256 3 DM particles and 256 3 gas particles in a box of 
comoving size L = 6.25 h^ 1 M pc. We perform s i mulat ions including photoion- 
isation by a uniform, evolving lHaardt Sz Madau ( 200ll ) UV background in the 
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optically thin limit at redshifts below the reheating redshiflQ z r . These simula- 
tions are denoted r[z T ]L6N256. To study the effect of photoheating, we compare 
these simulations to a simulation that does not include photoionisation and 
which is denoted L6N256. We explore a range of reheating redshifts to inves- 
tigate the sensitivity of our conclusions to changes in this parameter. We note 
that for the simulations that included the UV background the employed box 
size is sufficiently large and the employed resolution is sufficient high to ob tain 
converged results (for details see lPawlik. Schave. van Scherpenzeelll2009l ). 



3. The clumping factor 

Our main motivation for computing the clumping factor of the IGM is to evaluate 
the critical SFR density required to keep the IGM ionised. The critical SFR 
density describes a balance between the number of ionising photons escaping 
into the IGM (parametrised by / eS c) and the number of ionising photons that are 
removed from the IGM due to photoionisations of recombining hydrogen ions 
(parametrised by Cigm)- It is important to realize that only recombinations 
leading to the removal of ionising photons which escaped the ISM of the star- 
forming regions contribute to this balance. 

We employ a threshold d ensity criterion to separate the gas in t he ISM 
from the gas in the IGM fe.g.. lMiralda-Escude. Haehnelt. fc Reesll2000l h Ionis- 



ing photons are counted as escaped once they enter regions with gas densities 
Pb < Pthr- We treat the threshold density p t hr = Pb^thr as a parameter and 
compute the clumping factor C(< A t hr) = (p^) Ph < Pthv / 'p\ as a function of the 
corresponding threshold overdensity A tlir - This is done by performing a volume- 
weighted summation over all SPH particles with overdensities Aj < A t hr> i.e. 
we set C(< A thr ) = Ea^A^ W/£ a 

i<A thr hf, where hi is the radius of the 
SPH smoothing kernel associated with SPH particle i. 

By definition, C(< A t hr) increases monotonically with the threshold over- 
density Athr- We set an upper limit A* to A tlir , corresponding to the threshold 
density nj^ = 10 _1 cm~ 3 for the onset of star formation that we employ in our 
simulations. This choice is conservative, since the threshold density marking the 
escape of ionising photons and hence the clumping factor of the IGM to be used 
in Eq. [T]is likely lower. We conservatively identify the clumping factor C(< A*) 
with the clumping factor of the IGM, i.e. Cigm = C(< A*). 

The left-hand panel of Fig. [1] shows C(< Athr) f° r the simulations r9L6N256 
and L6N256 at redshifts shortly before (at z = 9.08) and well after (at z = 6) 
the reheating redshift z r = 9. The inclusion of photoheating strongly reduces 
the clumping factor for large threshold overdensities. This is because the asso- 
ciated increase in pressure smoothes small-scale density fluctuations. Note that 
because of the reduction of the clumping factor, photoheating provides a posi- 
tive feedback on reionisation, w hich acts in addition to the well-known negative 
feedbacks (for a review see, e.g.. lCiardi Fe"rr ara 2005). At z = 6 the clumping 
factor reaches a maximum value of C(< A*) ~ 6. This is five times smaller 



J If z r > 9, we use the z = 9lHaardt fc Ma dau fl200lT l UV background for all redshifts 9 < z < z r 
and employ the evolving lHaardt fc Madaul (2001) UV background for redshifts z < 9. 
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Figure 1. Left-hand panel: Clumping factor C(< Athr) of gas with over- 
densities A < A t hr- The inclusion of photoheating in r9L6N256 leads to a 
clumping factor that is substantially smaller than that obtained from L6N256, 
for threshold overdensities A t hr > 10. The maximum threshold overdensity 
we consider corresponds to the critical density njj = 10 _1 cm~ 3 (translating 
into an overdensity A*) for the onset of star-formation. Right-hand panel: 
Evolution of the clumping factor Cigm = C(< A*) for different reheating 
redshifts z r . At z — 6 we find Ctqm ~ 6, insensitive to the reheating redshift 
provided that z r ;> 9. Based on lPawlik. Schave. fe; van Scherpenzeell (|2009T ). 



than the value quoted in Gnedin &: Ostriker ( 19971 ). which has been employed 
in ma ny (observational) studies . The difference with our result is probably be- 
Gnedin fe Ostrikerl ( 19971 ) computed the clumping factor using a density 



cause 

threshold implicitly set by the maximum overdensity resolved in their simula- 
tion. The right-hand panel of Fig. [TJ shows the evolution of the clumping factor 
Cigm = C(< A*) for different reheating redshifts. It demonstrates that our 
estimate of the clumping factor of the IGM at z = 6, Cigm ~ 6, is insensitive to 
variations in the redshift at which the UV background is turned on. 



Discussion 



We computed the clumping factor Cigm = (Pb)iGM/ p\^ a measure for the aver- 
age recombination rate in the intergalactic medium (IGM), using cosmological 
smoothed particle hydrodynamics simulations. The clumping factor of the IGM 
depends critical ly on the definition of which gas is c onsid ered to be part of the 
IGM. Following iMiralda-Escude. Haehnelt. & Reesl ([2000), we assumed that all 



gas with densities below a threshold density constitutes the IGM and computed 
the clumping factor as a function of this threshold density. 

Even our most conservative estimate for the clumping factor, Cigm ~ 6, is 
five times smaller than the clumping factor that is usually employed to determine 
the capacity of star-forming galaxies to keep the z ~ 6 IGM ionised. Setting 
Cigm = 6 in Eq.[U the critical SFR density is p\ = 0.005 M yr" 1 Mpc~ 3 . 
This is smaller than recent observational estimates for the SFR density at z ~ 6, 
p\ ph 0.02 ± 0.004 M yr" 1 Mpc" 3 (integrated to the observed z pa 6 faint - 
end limit L > 0.04 L* z=3 and with new dust-correction; iBouwens et al.| [2009). 



provided that / esc > 0.25. Our study thus suggests that the observed population 
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of star-for ming galaxies may be cap able of keeping the z ~ 6 IGM ionised (see 
also, e.g., iBolton &; Haehneltll2007l ). In fact, recent studies suggest that, if we 
assume clumping factors as low as suggested by our work, observed galaxies may 
be capable of ke eping the Universe substantially ionised out to redshifts as high 
as z < 8 (e.g. 



Finkelstein et al. 1 120091 : lYan et al.l 120091: iBouwens et al.l [2009; 



Labbe et al. Il2009l : iBunker et al. lf2009l : iKistler et al.ll2009f )T 



The clumping factor i s an important ingred ient of (semi-)analytic treat- 
ments of reionisation(e.g., Barkana Lo eb 200l]). We have studied the im- 
pact of photoionisation heating on the clumping factor of the IGM assuming 
a uniform ionising UV background in the optically thin limit. In reality the 
photoheating process will, however, be more complex, with self-shielding be- 
ing only one example of the physical effects that we ignored. The validity of 
the approximations inherent to our simplified treatment should be assessed us- 
ing high-resolution radiation-hydrodynamical simulations of reionisation that 
include cooling by metals and molecules and also allow for feedback from star 
formation. For instance, kinetic supernova feedback has been shown to am- 
plify the_ejfecjs_ofphotoheating in evaporating the gas out of low-mass ha- 
los dPawlik Schavd 20091). which will af fect the clumping factor of the IGM 
(IPawlik. Schave. &: van Scherpenzeel I2009I L 
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